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Background  and  Objectives 


Background:  The  Navy  uses  hybrid  lidar- 
radar  for  underwater  detection,  ranging, 
communications,  and  imaging. 

-  Modulate  intensity  of  lidar  laser  signal  with  radar 
(RF)  waveforms 

-  Recover  radar  waveform  from  received  lidar 
optical  signal 

-  Apply  coherent  detection  and  other  radar 
techniques  to  process  received  signal. 

Objectives:  Enhance  hybrid  lidar-radar 
performance: 

-  Develop  and  evaluate  various  digital  signal 
processing  (DSP)  algorithms  that  will  enhance 
hybrid  lidar-radar  performance  in  a  variety  of 
underwater  environments. 

-  Implement  algorithms  via  DSP  hardware: 

•  dynamically  reconfigured  via  software  (accomplish 
multiple  missions  with  a  single  sensor) 

•  real-time  processing 

•  reduced  loss/temperature  sensitivity 
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Approach: 

•  Leverage  algorithms  that  are  widely  used  in 
radar  and  wireless  communications 

•  Spatial  filtering 

•  Frequency  domain  reflectometry 

•  Blind  signal  separation 

Challenges: 

-  The  underwater  propagation  channel  is  very 
different  than  the  above-water  environment 

•  Use  an  underwater  (UW)  optical  propagation  model, 
Rangefinder,  to  generate  simulated  data  to  test  DSP 
algorithms 

-  Many  COTS  DSP  hardware  platforms  (e.g. 
software  defined  radios)  are  suitable  for  above¬ 
water  radar  communications  but  lack 
performance  for  UW  light  detection  and  ranging 

•  Utilize  FPGAs  tailored  for  application 


Rangefinder 


Xilinx  Zynq  FPGA  Platform 


high  range  resolution, 


Performance  Goals 

Simultaneously  achieve: 
range  accuracy,  range  precision,  and  unambiguous  range 


1 .  Range  resolution 

-  Indicates  ability  to  resolve 
closely  spaced  targets 

2.  Range  accuracy 

-  Indicates  how  close 
measurements  are  to  the 
true  range  value 

3.  Range  precision 

-  Indicates  how  closely 
multiple  measurements  are 
clustered  together 

4.  Unambiguous  range 

-  Rmax  indicates  maximum 
distance  at  which  true  range 
is  determined  with  no 
ambiguity 
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Depiction  of  accuracy  and  precision 
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Low  Accuracy 
High  Precisions 

Desired  scenario: 

High  Accuracy 
High  Precision 


High  Accuracy 
Low  Precision. 


Low  Accuracy 
Low  Precision 
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Depiction  of  ambiguous  vs.  unambiguous  range 


Progress  and  Activity 


•  Project  start:  June  2011 

•  Year  1 

•  Participated  in  ONR  NREIP  program  at  NAVAIR 

•  Assisted  with  water  tank  experiments 

•  Characterized  software  defined  radios 

•  Identified,  validated  new  lidar  backscatter  reduction  technique:  spatial  filter 

•  Year  2 

•  Participated  in  ONR  NREIP  program  at  NAVAIR 

•  Spatial  filter  experiments 

•  Identified  new  ranging  technique:  frequency-domain  reflectometry 

•  Identified  new  backscatter  reduction  technique:  blind  signal  separation 

•  Summer  201 3  (planned) 

•  Participate  in  ONR  NREIP  program  at  NAVAIR 

•  Thorough  evaluation  of  new  techniques 

•  Validate  new  techniques  with  laboratory  experiments 

•  Publications  ! 

•  “Underwater  Laser  Rangefinder,”  Proceedings  of  SPIE  201 2,  Ocean  Sensing  and 
Monitoring  IV,  Volume  8372 

•  “Techniques  to  enhance  the  performance  of  hybrid  lidar-radar  ranging  systems,” 
Proceedings  of  MTS/IEEE  Oceans  2012 

•  “Technique  to  Extend  Unambiguous  Range  of  Hybrid  Lidar-Radar  Systems,” 
Proceedings  of  MTS/IEEE  Oceans  2013  (submitted) 
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Light  propagation  in  water 


Absorption  decreases  the  total  signal  level  at  the  receiver 


Scattering  degrades  image  contrast,  resolution,  and  reduces  range  accuracy 


Absorption  vs.  Scatter-Limited  Performance 


Absorption-limited  Scatter-limited 

detection  -  more  light,  detection  -  more  light, 

more  range  more ‘clutter' 


Although  absorption  and  scattering  are  two  separate  phenomena,  their  effects  on 
water  conditions  are  often  combined  together  into  a  single  parameter,  the 
attenuation  coefficient  c  (units:  m-1) 

Beam  attenuation  in  water  follows  an  exponential  decay  law,  where  cz  is  the 
number  of  attenuation  lengths  (a.I.) : 


P(c,  z)  =  P0  exp(-cz) 

After  traveling  one  attenuation  length,  the  optical  beam  is  attenuated  by  a  factor  of  1  -  ^  «  63.21% 


We  are  focusing  on  enhancing  performance  in  scatter-limited  scenarios 


Three  ways  to  improve  performance  in 
backscatter-limited  scenarios 

Simplified  lidar-radar  system 
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Radar  signal  generator 


1.  Wavelength  Selection  (A) 


Light  penetration 
open  ocean 


Light  penetration 
coastal  ocean 
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Radar  receiver 


2.  Modulation  Frequency  Selection  ( fm0d ) 
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Backscatter  Magnitude  vs  Modulation  Frequency 

Higher  mod olationj 
frequencies  are  better- 
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Modulation  Frequency  (MHz) 
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3.  Digital  Signal  Processing  (DSP) 


Our  focus  is  to  apply  DSP  algorithms 
to  the  receiver  platform  in  order  to 

improve  performance  by 
processing  the  received  signal. 
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Technical  Results 
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Goal:  automatically  detect  a  target  and  calculate  range 
without  any  a  priori  knowledge 


Calculated  range  (m)  Calculated  range  (m) 


Approach  #1:  Spatial  Filter  Results 
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Results  from  2012  Experiments,  c  =  2. 5  m 

Single-tone  ranging  with  fm=1 40  MHz 


-1 


4 

3.5 
3 

2.5 
2 

1. 


*  No  filter 
— Nominal  range 

- 

•  Spatial  filter 

1_ i_ i - 1 - 

i 

.5  2  2.5  3  3.5 

Nominal  range  (m) 

Dual-tone  ranging  with  ^=140  MHz,  f2=180  MHz 
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Simulation  of  Theoretical  Performance 

Spatial  filter  simulation  for  single-tone  at  140  MHz 
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•  Spatial  filter 
■  No  filter 
— Nominal  range 
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Nominal  range  (a.I.) 

Spatial  filter  simulation  for  dual-tone  at  140  and  180  MHz 
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Single-tone  140  MHz 


0.80  m 


2.65  cm 


;e  Improvement 
Spatial  Filtering 
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Dual-tone  140,  180  MHz 
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Goal:  automatically  detect  a  target  and  calculate  range 
without  any  a  priori  knowledge 


Ranging 


Backscatter  reduction 


Target  detection 
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Approaches: 

1 .  Single/Dual  Tone 


Spatial  Filter 


Phase  to  Range 


2.  Frequency  Domain 
Reflectometry 


Blind  Signal 
Separation 


Peak  Detection 


Frequency  (MHz) 


Approach  #2  Ranging  Algorithm: 
Frequency-Domain  Reflectometry  (FDR) 
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•  FDR  is  a  stepped  frequency  technique  developed  to 
detect  faults  in  RF  and  fiber  optic  cables 

•  Small  step  size  can  lead  to  large  unambiguous  range 
performance 

•  Requires  a  large  bandwidth  to  achieve  comparable 
precision  to  CW  methods 


■  ianalog  signal! 
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Sample  stepped  frequency  modulation 
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Parameter  Equation 


Unambiguous  range 


Range  precision 
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Target  to  backscatter  ratio  (dB) 


Approach  #2:  FDR  Simulated  Results 
with  Automated  Target  Detection 
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Simulation  of  FDR  ranging  performance  vs.  turbidity 
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Attenuation  lengths 

Simulation  of  FDR  performance  as  function  of  range 
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Simulation:  Fixed  target  range, 
varying  turbidity. 

•  BW  =1.2  GHz; 

•  Af=  1.17  MHz; 

•  8R  =  4.69  cm; 

•  Rmax  ~  48.0  m 

•  No  backscatter  reduction 
processing 

Results:  automated  target  detection 
works  to  13  a.I.  corresponding  to 
target-to-backscatter  ratios  of  -2  dB 

Next  step:  apply  blind  signal 
separation  to  reduce  backscatter 
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Amplitude  (a.u.) 


Approach  #2:  Blind  Signal  Separation 
for  Backscatter  Reduction 

Photons  reflected  from  objects  and  backscatter  in  the  scene  are 
recorded  as  a  single  "mixed"  lidar  return 
Applying  blind  signal  separation  (BSS),  the  target  and  scatter 
components  of  the  return  can  be  separated 

Frequency  response  without  BSS 
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Original  FDR  data  contains 
backscatter  and  target  signals 
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BSS  separates  the  signal  into  separate 
backscatter  and  target  signals 

Frequency  response  of  BSS  components 


Reduced 


Remove  undesired 
backscatter  component 
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Blind  Signal  Separation  (BSS) 
Math  Overview 


•  Objective:  Separate  observations  of  several  signals  into  the 
independent  signal  components  based  on  their  statistics 


•  Model:  X  =  AS 

X\  observation  matrix,  A:  "mixing"  matrix,  S :  source  matrix 

•  Methodology:  Select  a  statistical  property.  Adaptively  learn  a  weight 
matrix  W  ~  A -1  based  on  X.  W  should  maximize  separation  based 
on  the  desired  statistical  property.  At  convergence,  source 
components  are  estimated  as  S  =  WX  «  A~xX 


•  Challenges: 

•  Choice  of  statistical  property 

•  Choice  of  convergence  criteria 

•  Estimation  of  number  of  sources 

•  Choice  of  number  of  observations 

•  Computational  complexity 
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Approach  #2:  FDR  Simulated  Results  with 
and  without  Backscatter  (BSN)  Reduction 
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Target-limited 

Simulated  return  for  object  at  1  m 
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Target  >  BSN 

Simulated  return  for  object  at  3  m 
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Scatter-limited 


Simulated  return  for  object  at  7,5  m 
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Simulated  return  for  object  at  1  m 


Simulated  return  for  object  at  3  m 


Simulated  return  for  object  at  5.5  m 


Simulated  return  for  object  at  7.5  m 


Green  lines  indicate  actual  target  position 
Range  bin  size:  9.38  cm;  Turbidity:  c  =  1.6  m-1 


WlNlHt  Approach  #2:  FDR  with  BSS  and 

Automated  Signal  Detection 
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BSS  simulation  for  FDR  w  1 .2  GHz  BW  and  1024  tones 
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c  =  1.6  m_1;  BW  =  1 .2  GHz;  Af  =  1 .17  MHz 
SR  =  4.69  cm;  Rmax  =  48.0  m 


Real-Time  Implementation 


•  Several  basic  building  blocks  haye  been  implemented 
on  real-time  DSP  hardware: 


•  Baseband  demodulation 

•  Decomposes  RF  signal  into  two  orthogonal  vectors 

•  CORDIC  (Coordinate  Rotation  Digital  Computer) 

•  Given  two  vectors,  finds  the  magnitude  and  phase  angle 
of  the  resultant  vector 

•  The  phase  angle  can  be  used  to  calculate  range 

•  Instantaneous  frequency  demodulator 

•  Calculates  the  frequency  of  the  input  signal 


Example:  demodulation  of  sinusoidal  frequency  modulated  signal 

1 .  System  design  2.  Implementation  3.  Calculation 


Technique  Comparison 


Approach  1 

Single-  and  Dual-Tone  Ranging 

+  short  dwell  time 
+  very  high  precision 

-  single  range  result 

-  tradeoff  between  fmod  and  Rmax 

Spatial  Filter 
+  computationally  simple 
+  single  parameter 

-  logistically  complicated 

-  implementation  dependent  on  fmod 


Approach  2 

Freguencv-Domain  Reflectometrv 

+  large  unambiguous  range 
+  “range  spectra” 

-  large  bandwidth 

-  long  dwell  time 

Blind  Signal  Separation 

+  logistically  simple 
+  implementation  independent  of  fmod 

-  computationally  complicated 

-  many  parameters 
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Summary 


•  Significant  progress  has  been  made  in  developing  a  digital  signal 
processing  approach  to  enable  hybrid  lidar-radar  ranging  systems  to 
perform  automatic  target  detection: 

•  Ranging:  frequency-domain  reflectometry  technique  provides 

large  unambiguous  range  j 

•  Backscatter  reduction:  blind  signal  separation  technique  is  able 
to  separate  the  target  and  scatter  returns 

•  Target  detection:  simple  peak  detection  algorithm  is  able  to 
automatically  detect  target  position 

•  Planned  work 

•  Experiments  planned  for  summer  2013  to  test  performance  of 
new  technique  as  a  function  of  range  and  turbidity 

•  Continue  to  perform  experiments  with  spatial  filter  method 

•  MTS/IEEE  Oceans  Conference  presentation  in  Fall  2013 

•  Explore  more  sophisticated  target  detection  algorithms  that 
combine  FDR  and  signal-tone  approaches  to  improve  precision 

•  Develop  real-time  versions  of  each  technique 
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Frequency-Domain  Reflectometry 

(FDR)  Overview 


1 .  Transmit  stepped  frequency  $weep  with  N 
tones 

2.  For  each  frequency  ft,  record  magnitude  At  and 
phase  (pi 

3.  Construct  complex  vector  St  =  At  exp(/<Pi) 

4.  Calculate  inverse  Fourier  transform  of  Sf. 
s(t)  =  ^  /  S(f)  exp(2nft)  df 

5.  Peaks  appearing  in  the  time-domain  can  be 
converted  to  range  by  multiplying  by  the  speed 
of  light 
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Frequency-Domain  Reflectometry 
(FDR)  Range  Equation  Derivation 


FDR  range  equations  are  based  on  applying 
Nyquist  sampling  theory  in  the  frequency-domain 

Unambiguous  range: 

•  Frequency  sample  spacing:  A / 

•  Maximum  sampling  time:  tmax  = 

•  Convert  to  range:  Rmax  =  v-  tmax  =  ^ 

Range  precision: 


Time  bin  size:  St  =  —  L mny 

2 N  max  2NAf  2 B 

The  factor  of  one-half  is  needed  because  the 
maximum  error  is  half  the  size  of  a  bin 


•  Convert  to  range:  SR  =^St  — 
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Blind  Signal  Separation 
(BSS)  Math  Overview 


•  Objective:  Separate  observations  of  several  signals  into  the 
independent  signal  components  based  on  their  statistics 


•  Model:  X  =  AS 


X :  observation  matrix,  A:  "mixing" 


matrix,  S:  source  matrix 


•  To  separate  components  based  on  variance: 

1 .  Calculate  sample  covariance  matrix  C  =  XTX 

2.  Apply  singular  value  decomposition: ! 

X  =  USVT 

U\  projection  of  X  onto  eigenvectors  of  C,  S:  diagonal 
matrix  of  eigenvalues  of  C,  V:  eigenvectors  of  C 

3.  Remove  undesired  components  by  "zeroing-out"  columns  of 
U  and  corresponding  rows  of  S  to  obtain  U0  and  S0 

4.  Calculated  filtered  observations  as  XfM  =  U0S0VT 
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